Many-body band structure and Fermi surface of the Kondo lattice 
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We present a theory for the single particle excitations and Fermi surface of the Kondo lat- 
tice. Thereby we construct an effective Hamiltonian describing the creation and propagation of 
single-particle-like charge fluctuations on an 'RVB-background' of local singlets. The theory may 
be viewed as a Fermionic version of linear spin wave theory and is of comparable simplicity so that 
the calculations for the strong-coupling limit can be performed analytically. We calculate the single 
particle spectral function for the 'pure' Kondo lattice as well as for several extended versions: with 
a Coulomb repulsion between conduction and /-electrons, Coulomb repulsion between conduction 
electrons, and a 'breathing' /-orbital. In all cases we study the evolution of the spectrum in going 
from the Kondo insulator to the Heavy electron metal. We compare our results to exact diagonaliza- 
tion of small clusters and find remarkable agreement in nearly all cases studied. In the metallic case 
the /-electrons participate in the Fermi surface volume even when they are replaced by localized 
Kondo-spins and the number of bands, their dispersion and spectral character, and the nontrivial 
(i.e. non-rigid band-like) doping dependence including a pronounced transfer of spectral weight are 
reproduced at least semiquantitative^ by the theory. 
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I. INTRODUCTION 

The theoretical description of the Kondo lattice re- 
mains an outstanding problem of solid state physics. This 
model, or variations of it, may be viewed as the appropri- 
ate one for understanding such intensively investigated 
classes of materials as the heavy electron metals [Q^, 
Kondo insulators Q and possibly the recently discovered 

transition metal hydride-based switchable mirror com- 
poimds 1^,^ . The simplest model which incorporates the 
essential physics may be written as 

^ = E ^'^^..Cfc,. + + H.c.) 

k,a i.cr 
i,a i 



Here we consider the minimal model, where each unit cell 
contains two orbitals, one of them for the mobile con- 
duction electrons the other for the strongly correlated 
/-electrons. Then, cj^ ifjcr) creates a conduction elec- 
tron (/-electron) in cell i, and efc=-^ j e^'^'^^^~^^'> tij 
is the Fourier transform of the inter-cell hopping integral 
Uj for c electrons. For later reference we have already 
included a Coulomb repulsion U fc between / and conduc- 
tion electrons in the same cell, and a Coulomb repulsion 
Uc between conduction electrons. The latter two param- 
eters are usually taken to be zero, but as will be seen 
below our formalism allows to take them into without 
any additional effort. In the so-called symmetric case. 



Uf = 2ef and the limit U 
its strong coupling limit 



OO (^1), can be reduced to 



Hsc ^ E ^A,a^k,a + -^E -^^^^ ■ ^ 



(2) 



where Si^c i^ij) denotes the spin operator for conduc- 
tion electrons (/-electrons) in cell i and J = 4V'^/Uf. 
While the impurity versions of these models, which re- 
tain only a single /-site in a sea of conduction electrons, 
are well understood p|,p| ,p 12 and are even amenable 
to exact solutions |13 14[ |, much less is known about the 
lattice models. One problem which by many is believed 
to be at the heart of the solution is the way in which the 
more or less localized /-electrons, which in the strong 
coupling theory are replaced by mere spin degrees of 
freedom, participate in the formation of the Fermi sur- 
face and the heavy quasiparticle bands. Experiments on 
heavy Fermion compounds |l5| , computer simulations 
of Kondo lattices ||l^,|l^ and theoretical considerations 
|]l8| , p^ suggest that despite their frozen charge degrees of 
freedom, the /-electrons participate in the Fermi surface 
volume as if they were uncorrelated. In other words, 
the experimental Fermi surface volume corresponds to 
the case Uf = U/c = Uc = 0. The limiting cases V=Q 
or J—0, which obviously do not allow for participation 
of the / electrons in the Fermi surface, therefore repre- 
sent singular points, so that a perturbation expansion in 
the (small) parameters V or J may not be expected to 
give meaningful results. Rather, the interaction between 
/-spins and conduction electrons must be incorporated 
in a non-perturbative way, in a similar manner as the 
single-impurity Kondo effect |14|. It is the purpose of 
the present manuscript to present a minimum effort the- 
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ory for the Kondo lattice which is based on this basic 
requirement and shows how the nominal participation 
of the localized electrons in the Fermi surface can be 
understood even in the complete absence of any true 
hybridization. We describe the system by an effective 
Hamiltonian for the Fermion-like charge fluctuations on 
top of a strong coupling ground state, and show that this 
treatment leads to remarkable agreement with numerical 
results at least on energy scales which are relevant to 
high-energy spectroscopy. We would also like to point 
out that the method of calculation is similar in spirit to 
the cell-perturbation method developed by Jefferson and 
coworkers 12 0(] . A preliminary report has been published 
elsewhere [bill. 



II. CONSTRUCTION OF THE EFFECTIVE 
HAMILTONIAN 

As our starting point we choose the case of vanishing 
inter-cell hopping (i.e. ek = 0). The lattice prob- 
lem then reduces to single cell problems so that we first 
discuss the eigenstates of a single cell with 1, 2 and 3 
electrons. The two electron ground state is a singlet with 
wave function 



I* 



(2)\ 



) = [ afjfl + ^{c]fl + /|cl) + jc]c\ ] \vac) (3) 



/3 



(we have suppressed the site index on the Fermion cre- 
ation operators). The ground state wave function and 
energy is then obtained by diagonalizing the 3x3 matrix 

/~2ef + Uf , V2V , \ 

V2V , ~ef + Ufc , V2V . (4) 
V , V2V , Uc J 

For the strong coupling model (^) the problem becomes 
trivial with f3—l, the energy of the two-electron 

ground state is — ^ . 

The single or three electron states can be written as 

= ( + /3;'4cj/t ) \vac), (5) 

and the wave functions and energies are obtained by di- 
agonalizing the matrices 



Hi = 



-ef , V 
V , 



and 



'2ef + Uf + 2Ufc 



-V 

-£f + Uc + 2Ufc 



(6) 



(7) 



For the strong coupling limit the index v takes only one 
value and we have j—a— 0, and (3=P= 1. Both the single 



and three electron states have zero energy in this case. 
For later reference, we also define the photoemission 
(PES) and inverse photoemission (IPES) matrix elements 
(here a=c, f) 



,(2)\ 



(8) 



They can be expressed in terms of the wave functions 
defined above, e.g. 



aa, 



V2 



(9) 



where the additional sign in the first equation is due to 
our convention for ordering the two spin directions in (|^) 

H- 

We now return to the lattice problem and consider the 
case of half-filling (i.e. two electrons/unit cell, corre- 
sponding to the Kondo insulator). For vanishing ti^j the 
lattice ground state is simply the product of N single-cell 
ground states of the type (^) (see the state labeled (a) in 
Figure In the following, this state will be referred to 
as the vaccuum . Then, switching on tij produces charge 
fluctuations in the vacuum state: an electron can 
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FIG. 1. ' Chkrg^ fluctuations arid their' pr^ipagation (left 
panel) and their representation in terms of model Fermions 
(right panel). 

hop from cell i to another cell j, leaving the cell i in a 
single electron eigenstate with number u, and the cell j 
in a three electron eigenstate with number jj, (see state 
(b) in Figure |l|). In a further step an electron from the 
three-fold occupied cell j can hop to another neighbor j 
leaving cell j in a two-hole eigenstate and j in a three 
hole state (see Figure |l|c ) or, alternatively, an electron 
can hop from another neighbor i into cell «, leaving i in 
a two electron state, i in a single electron state. Finally, 
an electron from j can hop into i, leaving both i and j 
in two-electron states. In this picture the inter-cell ki- 
netic energy may be viewed of as a perturbation which 
has a two-fold effect: the pair creation of charge fluctu- 
ations and the propagation of these. It therefore plays a 
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completely analogous role as the transverse part of the 
Heisenberg exchange in the linear spin wave theory for 
the Heisenberg antiferromagnet, and in the following we 
want to exploit this analogy. To that end we make the 
additional assumption that a cell containing 2 electrons 

(2) 

must always be in the 2-electron ground state |5'g '). 
This means that after a charge fluctuation has passed 
through a given cell, the cell must return to the two- 
electron ground state. In other words, under this con- 
straint the propagating charge fluctuations do not leave 
a trace of excited cells - their propagation becomes coher- 
ent. Ways to relax this constraint and take into account 
incoherent processes will be discussed below. 
As a last important step, we note that the pair creation 
process in Figure ^ changes the energy by Ei, + E^ — 2Eq 
(which is nothing but the conductivity gap of a single 
cell), and that the switching of e.g. a hole- like charge 
fluctuation from species i' to species v' changes the en- 
ergy by e'^V — E^\ To keep track of these changes in 
energy, we interpret the difference fiy=E[P — e'^'^ as the 
'energy of formation' of the hole species v (and analo- 
gously for electron- like fluctuations). 
We now define our restricted set of basis states: 



\niVi,n2V2, 



■nNVN) 



1=1 



(10) 



with the side condition that = 2 automatically im- 
plies Vi = 0. In the following, we will diagonalize the 
Hamiltonian in the subspace of the states (|l^) . To that 
end, we represent the basis states ([l0| ) in terms of 'model 
Fermions' j2^: if a cell with number i is in the two- hole 
ground state we say it is empty; if the cell is in the i^*'* 
single electron state with z-spin a we model this by the 
presence of a hole-like Fermion, created by a\ ^; and 
if there are three electrons forming the /i*'' single cell 
state with z-spin a we say that the cell is occupied by 
an electron-like model Fermion, created by b\ ^ ^. Then, 
solving H in the restricted basis (|l^) obviously is equiva- 
lent to diagonalizing the Hamiltonian Rgj f=VHV where 



•a ,o" 



(11) 



with 
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Here V projects onto the subspace of states where no site 
is occupied by more than one Fermion. This kinematic 
constraint reflects the fact that the state of a given cell 
must be unique. Due to the product nature of the basis 
states (|l^) , the evaluation of the matrix elements of the 
intcr-ccU kinetic energy (^ reduces to the calculation of 
matrix elements between products of no more than two 
single cell states. The matrix elements on the r.h.s. of 
( p^ ) are therefore simply products of the c-like photoe- 
mission and inverse photoemission matrix elements for 
a single cell (^; this corresponds to the quite intuitive 
picture that the propagation of a c-electron is equivalent 
to photoemisson in one cell, and inverse photoemission 
in the neighboring one. The respective matrix elements 
give the renormalization of the inter-cell hopping due to 
intra-cell (i.e. local) correlation effects. Also, as long as 
the interaction between electrons contains only intra-cell 
terms the entire strong correlation physics obviously is 
completely taken care of by the calculation of the single- 
cell states, and only enters via the single cell energies 



E^^ and 



E, 



(3) 



E, 



(2) 



Thus, while we 



are presently only using Coulomb repulsions as intra-cell 
interactions, Hund's rule exchange or electron phonon 
coupling could also be treated in the same way. 
Having computed the matrix elements and excitation en- 
ergies the most obvious next step then is (in analogy 
to linear spin wave theory) to relax the constraint en- 
forced by 7^, whereupon the Hamiltonian ( ITTl ) is readily 
solved by Bogoliubov transformation. This gives us the 
energies and dispersion, and for the full Kondo lattice 
Hamiltonian we obtain 4 bands (we have two hole-like 
and two electron- like model Fermions), whereas for the 
strong coupling version we have only two. For the latter 
case, (0) takes the form 



- J E ''9n{cT) tk (4,.alfc,<. + H.c). (13) 



This is readily solved by the ansatz 



(14) 



and, introducing A=3J/2, we obtain the quasiparticle 
dispersion 



i?±(fc) = (l/2)[efc±./6| + A 



(15) 



shown in Figure At half-filling, particle-hole symme- 
try requires the chemical potential to be zero, so that 
the lower of the two bands (|l^) is completely filled, the 
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upper one completely empty. We note that formally ( |l5| ) 
is completely equivalent to the hybridization of a disper- 
sionless effective /-level in the band center with a free 
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FIG. 2. Quasiparticle dispersion (top panels) and disper- 
sion of the c-like spectral weight along the lower band (bot- 
tom panels) for the ID strong coupling model with J = 0.5, 
Efe = — 2cos(fc). The density of conduction electrons is f in 
(a) and 0.8 in (b). 

electron band with dispersion e^, the strength of the 
nominal mixing element being A. This results in the 
heavy , i.e. almost dispersionless bands immediately 
above and below the Fermi energy in Figure ||- the slope 
of these bands decreases strongly with decreasing J. It 
should be noted, however, that the resulting energy gap 
of A does not arise from the formation of a bonding and 
antibonding combination of c-like and /-like Bloch states, 
as in the hybridization model; rather, this gap originates 
from the energy cost to break two intra-cell singlets in 
the first step of a charge fluctuation. This gap therefore 
is of a very similar nature as the energy gap in a super- 
conductor: the minimum energy for moving an electron 
from some site i to a distant site j (which by definition 
is the single particle gap of the system) is two times the 
energy required to break a pair (which may be either a 
Kondo singlet or a Cooper pair: the flrst pair is broken 
at site z, because one member of the pair is removed - 
the second pair is broken at site j because the surpuls 
electron interferes with the pair formation around this 
site. 



III. SINGLE PARTICLE GREEN S FUNCTION 

To compute the full single particle spectral function we 
need to resolve the ordinary electron creation and annihi- 
lation operators in terms of the model Fermions, a and h. 
Taking into account our basic assumption, namely that 
a single cell with 2 electrons can only be in its ground 
state, we can expand the electron annihilation operator 
(where a = c, f) 



Taking matrix elements of both sides we readily find 
that Uv=ra,v,a, v^=s*^^^^^. Next, we can replace e.g. 

\^^^\) (^'o^^ I ^ a| y and thus have the desired resolution 
of the photoemission operator. Specializing to the strong 
coupling limit, the annihilation operator for c-electrons 
takes the form 



(17) 



Using ( p^ we can now resolve the annihilation operator 
in terms of the quasiparticles, obtain its matrix element 
and square it to obtain the spectral weight in the lower 
(i.e. occupied) band as 



W 



1 



(1 - 2sign{(j)uk,aVk,a) 



(18) 



Let us assume that ^ A and let k be in the outer part 
of the Brillouin zone, i.e. we assume that we are deep in 
the heavy portion of the occupied band (see Figure |^) . 
Then, we find 

(A/6fc)2 « 1 

i.e. the heavy band does have an extremely small spec- 
tral weight (see the lower part of Figure |^a) . This will 
turn out to be of considerable importance in a minute. 
To proceed to the doped case, we need an expression for 
the electron number operator. While at first sight this 
may appear a triviality, we will now see that one thereby 
runs into a rather deep-rooted problem, which reflects 
the special features of the strong correlation problem. In 
the vacuum state the number of electrons (counting c 
and /-electrons) is 2N and the presence of an a-Fermion 
(6-Fermion) decreases (increases) the electron number by 
1, so that for the full Kondo lattice the electron number 
operator should be simply 



E E 7L,^7fc,M,^ - 27V. 



(19) 



For the strong coupling limit, we obtain in an analogous 
fashion 



(20) 



The extra —2N on the r.h.s. of ( p9| ) simply cancels the 
lowest of the 4 bands obtained in the full Kondo lattice. 
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As will be seen below, this is a practically dispersionless 
lower Hubbard band for the /-electrons, which is absent 
in the strong coupling limit (or better: pushed to -co). 
The Fermi surfaces of both models thus are completely 
equivalent, and e.g. for the case of hole doping in ID one 
would obtain a Fermi momentum kF=^iPc + 1) (where 
Pc < 1 denotes the density of conduction electrons) in 
the lower hybridization band. This implies that we have 
a Fermi surface which satisfies a nominal Luttinger the- 
orem, i.e. the / -electrons are treated as participating in 
the Fermi surface. The physical origin, however, is the 
fact that we have a density 1 — pc oi holes in the singlet 
background , which forms the vacuum for our treatment. 
This results in a hole pocket centered on fc = tt with 
Fermi momentum f (1 — Pc)- The resulting Fermi surface 
then is nominally equivalent to a Luttinger Fermi surface 
including the /-electrons and obviously this equivalence 
holds true irrespectively of dimensionality and Fermi sur- 
face topology. We believe that this is the reason why 
the /-electrons seem to participate in the Fermi surface 
volume despite the fact that they are localized. At half- 
filling, we do not have a half-filled band of single-particle 
like mixtures of c and /-electrons - such a picture is ob- 
viously completely wrong for the strong-coupling model. 
Rather, the half-filled ground state should be viewed as 
an array of local singlets, where each /-spin captures one 
conduction electron to form an immobile singlet. The 
system has a single particle gap because removing an 
electron and re-inserting it at a far away site breaks two 
of the local singlets, in a completely analogous fashion 
as in a superconductor. Unlike a superconductor, the 
two members of each pair do belong to different species , 
whereby the number of one of the species, the /-electrons 
is essentially independent of the total electron density. 
Removing electrons from the system thus produces holes 
in the singlet background and leaves behind unpaired /- 
electrons. Unlike the actual /-electrons, which are largely 
immobile (completely immobile for the strong-coupling 
model) the holes are mobile, because an /-electron can 
form a pair with a c-electron from a nearby pair, thereby 
leaving another /-electron unpaired. 
However, there is still a major complication: the electron 
number must also equal the k and lu integrated photoe- 
mission weight, which is given by the expectation value 
of the operator 



-sign{<7){bl^^a^ 



k,(j) 



(21) 



Inserting the expansion ([T^) of the c and / electrons in 
terms of the model Fermions into this expression it is 
easy to see that in general N^. For example in the 

strong coupling limit we find using ( p!?] ) 



k.ij"--k,S 



-krM]+N. (22) 



We thus arrive at the at first sight devastating conclusion 
that counting the electrons in real space on one hand and 
integrating the spectral weight in fc-space on the other 
hand give us different results for the electron number. 
On the other hand, this is not an indication for a quali- 
tative flaw in our theory, but has a very clear and simple 
physical origin, namely the fact that in a strongly corre- 
lated electron system spectral weight and band structure 
are completely decoupled. As an example, let us consider 
the band structure shown in Figure |^a. If one simply 
were to introduce a chemical potential corresponding to 



the real-space electron count (20) the Fermi energy would 
necessarily cut into the heavy part of the lower band. Let 
us assume that the only effect of doping were that kp pro- 
gressively cuts deeper into the heavy band portion of the 
otherwise completely rigid quasiparticle band structure 
calculated above. Trivially, for the II? chain the number 
of occupied momenta, A^occ would change with kp ss, 



ON, 



dki 



= 4. 



On the other hand, the integrated photoemission weight, 
NpES, would change as 

okp ekp 

Introducing a simple chemical potential following (|2C 
into an otherwise rigid quasiparticle band structure thus 
inevitably leads to the breakdown of the sum rule for 
the integrated photoemission weight. The breakdown of 
rigid-band behaviour and the failure of a simple electron 
count are therefore a completely natural consequence 
of the special feature of strongly correlated electrons, 
namely that they may have quasiparticle weights which 
substantially deviate from unity. 

In order to cope with this problem, we try the simplest 
possible solution and enforce the consistency of ordinary 
electron count and spectral weight integration by adding 
both expressions, (po| ) and (|2^), for the electron number 
to the Hamiltonian, each one with a separate Lagrangian 
multiplier: 



H ^ H - /iiVe - \N'^. 



(23) 



The notion of two chemical potentials may seem awkward 
at first sight, but as will be shown now, this approach 
results in a remarkable consistency with the numerical 
results. Let us again consider the strong coupling limit. 
The spectral weight operator Nf, takes the same form as 
the kinetic energy, but with the replacement — > 1 or, 



equivalently, U 



1. The Hamiltonian thus becomes 



1 



fc,£T 



k.CT 



3J 

T 



M 4 



k,a"'k,a 
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3J 



It.'' 



gn{a){ek - A) (5^ 



H.C.), (24) 



k.{T 



SO that, using again (p^), wc find the dispersion 

E±{k) = l[(ek~X)± - A)2 + A2 ] - pi. (25) 

Using the representation of the spectral operators, we 
obtain the momentum distribution/spin direction of the 
conduction electrons: 



c 1/1 efc-A 



2^ ^(efc-A)2 + A2- 



(26) 



Then, for |efc — A| 3> A we may replace the denom- 
inator on the r.h.s. by —\ek — A| and obtain — 
i(l — sign{ek — A)). This is simply the free-electron re- 
sult. For efe — A = on the other hand, we find n% = 1/2. 
Let us now discuss these results. To begin with, the real 
space chemical potential fi acts like a standard chemical 
potential, which cuts into the heavy band and, as dis- 
cussed above, produces a Fermi surface consistent with 
the nominal Luttinger theorem. On the other hand, the 
chemical potential for the spectral weight. A, gives rise 
to a pseudo Fermi surface for the conduction electrons, 
where drops sharply but continuously from a value 
« 1 to nearly (see Figure ^a). Since the integrated 
c-distribution must equal the number of conduction elec- 
trons, it is clear that A « ^{kp), with kp the Fermi mo- 
mentum for unhybridized conduction electrons. Carrying 
on the formal analogy of ( p5|) with a hybridization gap 
picture, A obviously plays the role of an on-site energy 
of the effective / level, and therefore the line = X 
marks the locus in k space where the strongly dispersive 
conduction-band bends over into the nearly fiat heavy 
band (see Figure |^) In other words, the band structure 
of the Kondo lattice is equivalent to an effective /-level, 
which is pinned near the frozen core Fermi energy for con- 
duction electrons of density pc, and which mixes into the 
conduction band with a matrix element of strength 3 J/4. 
We note that this is very much what one would expect 
intuitively: in the limit i 3> J the kinetic energy of the 
c-electrons is by far the dominant energy contribution of 
the system. This can be expressed as 



{Hkin) = J dknl efc 



and obviously the minimum value compatible with the 
Pauli principle is obtained by the free electron distri- 
bution for the conduction electrons. Then, in order to 
recover a (small) additional energy ~ J from the Kondo- 
hybridization, the system will not sacrifice much kinetic 
energy, i.e. n^(fc) will stay close to its free-electron shape. 
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FIG. 3. Momentum distribution for the strong coupling 
limit for different values of J, pc = 0.5. 

Figure ^ shows nc(fc) for a ID chain of the strong cou- 
pling model with pc—0.5, evaluated numerically by solv- 
ing self-consistently for p and A, for different values of 
J/t. For very large J (which is an unphysical limit- 
ing case) nc{k) is nearly constant and drops to at 
kp — 37r/4, the Fermi momentum for hybridized conduc- 
tion and / electrons. As J is reduced, the Fermi surface 
discontinuity shrinks more and more (but stays finite at 
any J) and the pseudo Fermi surface at /c^=7r/4 starts to 
develop. This behaviour of nc{k) is in almost quantita- 
tive agreement with the DMRG results of Moukuri and 
Caron |ll6|] , which may show the quality of our simple 
analytical calculation. 

Having discussed the strong coupling version, we now 
turn to the full Kondo lattice problem. Thereby we en- 
counter a new problem. In the strong coupling version, 
we have introduced the pseudo chemical potential A in or- 
der to enforce the consistency of real-space electron count 
and integrated photoemission weight; clearly we will have 
to do the same thing for the full Kondo lattice. For the 
full model, however, there is an additional distinction to 
be made because now we have c-like and /-like spectral 
weight, and we need to make sure that the total spec- 
tral weight is distributed between the two species in a 
proper version. We will see that this enforces to intro- 
duce yet another Lagrange multiplier which essentially 
governs the ratio of c-like to /-like weight. 
To that end, let us first construct an operator which 
counts the total number of a-type electrons (where a = 
c, /). To begin with, we define the expectation values of 
the electron numbers in the single cell states: 

These expectation values are readily computed from the 
single cell wave functions, e.g. n^^'^''=2a'^ + 0^ . By anal- 
ogy with ( |2^ ) we may thus write down the following op- 
erator to count the total number of a-electrons: 

+ (27) 



6 



Usin£ 



4*'"^ +n^''^=i we find Nc+Nf=Ne, with given 
by (19) (as it has to be). On the other hand by using 
the resolutions (|l6|) of the a operators we can form the 
operators N'^=J2k a ^fc a'^^.a, which give the integrated 
spectral weight for the a-electrons alone. We then re- 
place the Hamiltonian 



H 



H-f,N,~XiN'^ + N'f)-XfiNf-N'f). (28) 



We determine A from the condition {N^ + Nj.-)~Ne (A 
thus has precisely the same meaning as for the strong 
coupling limit discussed above) and A / from the require- 
ment {Nf — N'j:) — 0. This implies that automatically 
also {Nc — N^) = 0, whence we have reached 'species 
wise' consistency between real-space electron count and 
integrated spectral weight. 
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IV. COMPARISON WITH EXACT 
DIAGONALIZATION 



Using the eigenenergies obtained by diagonalizing the 
effective Hamiltonian (^Sj) with the self-consistently de- 
termined values of the Lagrange multipliers A and A/ 
we obtain the band dispersions. Combining the resolu- 
tions ( p^ ) of the c and /-operators in terms of the model 
Fermions and the eigenvectors obtained by diagonalizing 
( p8| ) we can compute the photoemission and inverse pho- 
toemission matrix elements, so that we can obtain the 
photoemission spectrum 



1 



U! - 



{H~Eo 



^0^ 



and the inverse photoemisison spectrum 



'io-iH-Eo) 



where a — c,f. Thereby |1'o) denotes the ground state 
wave function, Eq the ground state energy. We have 
carried out this calculation for several 1-dimensional ver- 
sions of the Kondo lattice and in the following we com- 
pare them to the results of Lanczos diagonalization. For 
the Lanczos studies we used a 6-unit cell chain. To sim- 
ulate longer chains, we combine spectra calculated with 
periodic and antiperiodic boundary conditions (see Ref. 
I prf for a detailed discussion) . While there is no rigorous 
justification for this procedure, inspection of the numeri- 
cal spectra shows that one can obtain remarkably smooth 
band structures in this way. 

To begin with. Figure ^ shows the Lanczos result for the 
single particle spectral function for a ID chain of the 
'pure' Kondo lattice, i.e. Ufc=Uc= 0, Figure || shows the 



FIG. 4. Single particle spectral function for the ID Kondo 
lattice, obtained by Lanczos diagonalization of a 6-unit-cell 
system. Full lines (dashed lines) correspond to c-like (/-like) 
spectral weight. The vertical dashed line gives the Fermi en- 
ergy Ef, peaks to the right (left) of this line correspond to 
electron creation (annihilation). Parameter values are f/ = 8, 
ej =4,V = 1, e{k) = -2cos(fc). 



result obtained from the theory (similar Lanczos spec- 
tra with slightly different parameters have also been ob- 
tained by Tsutsui et al. JitI). To begin with, in contrast 
to any band theory approach, our theoretical spectrum 
correctly reproduces the 4 well-distinguishable bands in 
the numerical spectra: the practically dispersionless up- 
per and lower Hubbard band , which have almost pure 
/-character, and the two hybridization bands, which 
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FIG. 5. Theoretical spectrum for the same parameter val- 
ues as Figure 0. 
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resemble the strong coupling- result The basic idea 

of our approach, namely to broaden the ionization and 
affinity states of a single cell into bands thus obviously 
works quite well. The spectra nicely show the Pseudo 
Fermi surface for the c-electrons: while the true Fermi 
momentum in the doped case is at fci?=57r/6, where the 
heavy band intersects Ep (although this is not easily rec- 
ognized in the theoretical spectra), there is a pronounced 
transfer of c-like spectral weight from below to above Ep 
at fc^ = tt/2. This means that doping shifts the drop of 
the momentum distribution for the c-electrons (which is 
nothing but the integrated c- weight below Ep) from 
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tt/2 to 7r/3 - this is preciseley the pseudo Fermi surface 
discussed above for the strong coupling limit (see Figure 
H). The theoretical spectra somewhat overestimate the 
weight of the heavy /-like portions in the hybridization 
bands, moreover the position of the /-like Hubbard bands 
is shifted to slightly too high (binding) energies. On the 
other hand, the Fermiology is reproduced quite well. It 
is interesting to note that qualitatively identical results 
for smaller values oi Uf and e / and lower electron filling 
have been obtained by Tahvildar-Zadeh at al. ||2^ using 
the D oo technique - quite obviously there is overall 
agreement between different numerical approaches, for 
quite different parameters values and fillings (and our 
analytical results). 

Figure ^ shows a more asymmetric case, more precisely 
we set £/=[/• As expected, the upper Hubbard band for 
the /-electrons is virtually absent. Interestingly, the the- 
oretical spectrum shows a very weak band at ~ 3t above 
the Fermi energy, whose intensity shows a weak increase 
with doping. The Lanczos spectrum for half-filling shows 
some diffuse /-weight roughly in this area, and even some 
indication of peaks in the doped case. Next, the disper- 
sion of the hybridization bands shows a quite pronounced 
asymmetry between photoemisison and inverse photoe- 
mission, which is nicely reproduced by the theory. The 
explanation is simple: due to the reduced energy of the 
pure /-state /|/||wac) the weight of this state in the 

single cell ground state for two electrons I^'q ) must in- 
crease. The weight of the two remaining states there- 
fore must decrease, so that the overall probability to find 
a c-electron in the ground state is reduced (this man- 
ifests itself also by the significantly smaller integrated 
c- weight in the photoemission spectrum). Since only the 
c-electrons can hop between cells, the band-width seen in 
photoemission therefore will be reduced. Next, since the 
c-weight which disappears from the photoemission part 
must reappear in the inverse photoemission spectrum, it 
follows from analogous considerations that the dispersion 
in the IPES part will be enhanced. As was the case in 
the symmetric case, doping causes a shift of the chemical 
potential into the lower hybridization band, but there is 
now also some change of the spectral character of the up- 
per hybridization band near k — tt: there the c-character 
increases, the /-character decreases with doping, and the 
theoretical spectra obviously do correctly reproduce this 
trend. 

Next, we study a more complicated form of the hopping 
term for conduction electrons. More precisley, we intro- 
duce a hopping integral between second nearest neigh- 
bors, and to have a pronounced effect we choose this to 
be of equal magnitude, but opposite sign as the nearest 
neighbor hopping. In other words, the conduction elec- 
trons now have the dispersion relation 



FIG. 7. Theoretical spectrum for the parameter values as 
in Figure H. 
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-2i(cos(fc) - cos(2fc)). 



(29) 



8 




3 -4 -2 





2 4 e 


1 1 1 1 1 1 
<n>=5/3 


I ' 

1 

,1 , 


1 1 1 1 1 1 1 1 

L A 




.J^L/^.L_ A...... 




l: A 
















, 1 1 , , , 1 1 , , , 1 , y.-;.V , 


-4 -2 





2 4 6 



FIG. 8. Lanczos spectrum for conduction electrons with 
the dispersion , all other parameter values as in Figure m 



For the band fillings under consideration, a system of 
noninteracting electrons thus would have 4 Fermi points 
rather than 2, and we want to see if our simple 'rule 
of thumb' deduced above for the band structure in the 
strong coupling limit continues to be valid also in this 
more complicated situation. Namely the band structure 
of the Kondo lattice should be qualitatively given by the 
one obtained for mixing with a dispersionless /-level, that 
is pinned to the 'frozen core' Fermi energy of the non- 
hybridizing conduction band, see Figure ^. Then, figure 
^ shows the Lanczos result for the single particle spectral 
function and Figure ^ shows the theoretical 
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FIG. 9. Theoretical spectrum for the parameter values as 
in Figure H. 



data. It is quite obvious that the theory remains valid 
also in this case. The relatively complicated band struc- 
ture predicted by the rule of thumb indeed can be seen 
quite well in the Lanczos spectra: the s-shaped conduc- 
tion band, with a local maximum at fc = 0, which crosses 
from inverse photoemission and back again and thereby 
mixes with the dispersionless, low intensity /-band (it 
should be noted that the low intensity /-like peaks at 
5n/6 and tt in the Lanczos spectra for (n) = 2 are actually 
on the photoemission side) . With doping away from half- 
filling the heavy band near fc = tt develops a Fermi edge, 
and the two Pseudo Fermi surfaces for the c-electrons at 
7r/6 and 27r/3 shift as expected for free, nonhybridizing c- 
electrons. In the doped case, the fiat low intensity band 
which skims below Ep and crosses Ep at kp — bir/G 
can be identified particularly well. Also the heavy part 
in the inverse photoemission at k — tt/3 and tt/2 which 
is predicted by the theory can be identified reasonably 
well. Again, the theoretical spectra somewhat overesti- 
mate the spectral weight of these heavy band portions, 
and the relative position of the lower Hubbard band is not 
correct. Our simple picture of the Fermiology, however, 
obviously remains valid also with this more complicated 
Fermi surface topology. 

We proceed to the case of non-vanishing Ufc, i.e. a 
Coulomb repulsion between / and conduction electrons. 
The Lanczos spectra for this case are shown in Figure 
the theoretical spectra in Figure The main effects of 
Ufc is to open a wider gap at half-filling, to enhance 
the weight of the heavy bands and to shift the Hubbard- 
bands to higher energy. All in all, the total width of the 
spectrum is now 12 = Uf + Ufc rather than Uf as it 
used to be in the preceding cases. As compared to the 
preceding cases, the Hubbard bands moreover acquire an 
appreciable dispersion, and doping causes a pronounced 
spectral weight transfer from the upper to the lower Hub- 
bard band. The widening of the gap at half-filling ap- 
pears hard to understand at first sight, because Ufc by 
itself does not increase the energy cost for a charge fluc- 
tuation: transferring a c-electron from one cell to another 
(see the process Figure |l|a ^d) does not change the en- 
ergy due to f-c repulsion as long as the /-occupation is 
nearly constant at 1 (which is almost certainly the case 
for the large we are using). It is straightforward to see, 
however, that the opening of the gap is due to a loss of ki- 
netic energy: switching on Ufc increases the energy of the 
state --i=(c|/| + /|c|) by Uf^ so that the energy differ- 
ence relative to the other two basis states is decreased by 
this amount (see the Hamilton matrix (|^)). This means 
that charge fluctuations in the two-electron ground state 
are enhanced, whence the kinetic energy for this state 
becomes more negative. On the other hand, this mech- 
anism for enhancing the charge fluctuations is operative 
neither for one nor for three electrons in a cell (see the 
respective Hamilton matrices (pl|) 
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FIG. 10. Lanczos spectrum with U /c=4, all other parame- 
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FIG. 12. Lanczos spectrum with I7c=4, all other parameter 
values as in Figure ^. 



and , so that there is no gain in kinetic energy in these 
cases. The charge fluctuation thus is accompanied by a 
net loss of kinetic energy, which in turn results in a larger 
charge gap. Moreover, it follows from the expression ( |2^ ) 
for the c-type spectral weight of the lower hybridization 
band that the distance in fc-space over which the weight 
drops is approximately given by Ac/vp, with the 
charge gap and vp the Fermi velocity. It is thus im- 
mediately obvious that an increase of the charge gap (for 
whatever reason) will give a 'more homogeneous' c-weight 
along the hybridization bands. While at half- filling the 
agreement between Lanczos and theoretical 
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FIG. 11. Theoretical spectrum for the parameter values as 
in Figure [iG. 



spectra is very good, the situation changes for the doped 
case. The lower of the two hybridization bands is still 
well described, but the weights and dispersions of the 
other bands, as well as the weight shift between the 
Hubbard bands are not reproduced well. Actually the 
Lanczos spectra show a rather strong reduction of the 
gap between the two central bands upon doping - we 
believe that this indicates a rather profound doping- 
induced reconstruction of the electronic structure and 
since our theory corresponds more or less to an 'expan- 
sion around the Kondo insulator' it may not be expected 
to reproduce this. 
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The situation is actually quite similar if we introduce a 
Coulomb repulsion between conduction electrons [p5|-p7| , 
Uc, see Figure 12 and |l^. Again, this opens a wider 
quasiparticle gap in the spectrum and the hybridization 
bands now have a predominant c-character over their en- 
tire width. At half-filling the upper hybridization band 
intersects and mixes with the upper Hubbard band for 
the /-electrons. The upper and lower Hubbard band 
for the /-electrons remain unaffected and are dispersion- 
less, the lower one has practically pure /-like character. 
At half-filling, there is again good agreement between 
Lanczos and theory but for the doped case the situation 
is again different. The Lanczos spectra show a rather 
dramatic collapse of the gap between the hybridization 
bands, the upper hybridization band, which intersected 
the upper Hubbard band at half-filling is now at least 
2t below. Again, this suggests a strong doping induced 
reconstruction of the entire electronic structure, and our 
theory naturally fails to reproduce this. 
The effect of Uc at half-filling can be understood already 
in the context of the strong coupling model limit, where 
we simply would have to replace A — > A -I- f/c- The more 
homogeneous c-character of the hybridization bands then 
follows in an analogous fashion as for Ufc- 
As a last example, we turn to a somewhat exotic ver- 
sion of the Kondo lattice model, where the magnitude of 
the c-f hybridization depends on the occupation of the 
/-level. More precisely, we replace 

with Vi ^ V2. Such a model may be relevant |Q to 
describe the recently discovered Q transition metal hy- 
drides with switchable mirror properties, such as YH3. In 
this case, the Yttrium 4(i electrons would play the role of 
the conduction electrons , the Hydrogen corresponds to 
the /-electrons . The conditional hopping term is sup- 
posed to describe the relaxation of the orbital wave func- 
tion on Hydrogen as a function of electron occupation 
0. This is manifested e.g. by the dramatically larger 
radius of the free ~ ion as compared to the neutral free 
H atom (for a detailed discussion see Ref. ||^, for a dis- 
cussion of such an 'orbital Kondo effect' in the context of 
cuprate superconductors see Refs. p^ , p9t ). Here we are 
not so much interested in the details of the correspon- 
dence with YH3, but rather in the applicability of our 
theory to this model. The conditional hopping is easily 
incorporated by replacing the Hamiltonian for two 
electrons in a cell by 



Ho 




'£/ + Ufc 




(30) 



and analogous replacements in the single and three- 
electron subspaces. It is quite obvious, that the con- 
ditional hopping will lead to a dramatic increase of the 




FIG. 14. Lanczos spectrum for the breathing / orbital . 
Parameter values are U — ef — 2, Vi = 2.0, V2 = 0.2, 
ek = — 2cos(fc). 

charge gap: in our picture, the magnitude of the gap is 
determined by the energy difference between two cells 
with 2 electrons on one hand, and a cell with 1 and one 
with 3 electrons on the other hand. Then, the electron 
in the singly occupied cell can mix with the conduction 
band only by using the (much smaller) hybridization in- 
tegral V2, which corresponds to the collapsed / orbital. 
Only the single hole in the cell occupied by 3 electrons 
can delocalize using the large hopping integral Vi. A 
charge fluctuation thus will result in a huge loss of ki- 
netic energy and thus open a substantial gap in the band 
structure 
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(the effect is thus similar, but much stronger than the 
Ufc discussed above). This gap will be opened even in 
the complete absence of any Coulomb repulsion on the /- 
orbital, so that the breathing /-orbital can lead to strong 
correlation-like behaviour even if there is no really strong 
Coulomb repulsion at play. In fact, comparing the Lanc- 
zos and theoretial results (Figures O and 15) one can 



note first of all a good agreement between the two, and 
in addition a very substantial correlation gap in the spec- 
trum. Also, there are low-intensity 'Hubbard bands', 
whose position and spectral weight are reproduced quite 
well by the theory. An interesting feature in the Lanczos 
data is the extra /-like band in the inverse photoeniission 
spectrum which disperses from just above Ep at fc=7r to 
sa 2t above Ep at fc=0, and which is completely absent in 
the theory. We believe that this band is a consequence of 
superconducting pairing correlations in the doped ground 
state. As discussed by various workers, the breathing /- 
orbital presents a very strong pairing mechanism [ p8|j29| 
for electrons in the doped material, in that the electrons 
will always pair up around one /-site to hybridize by the 
large hopping parameter Vi . In the spectral function the 
presence of such pairing correlations, will manifest itself 
as a Bogoliubov-type 'mirror image' of the topmost band, 
i.e. the replacement 



fife - M 



where Ck is the 'nonsuperconducting dispersion' of the 
lower hybridization band (which should correspond to 
the one given by our theory) and Age the superconduct- 
ing gap. The latter appears to be quite small in the 
numerical spectra. We defer a detailed discussion of su- 
perconductivity to a separate publication, but wc note 
that our theory does not encompass superconductivity, 
so that the absence of this band in the theroretical spec- 
tra can be no surprise. 

Summarizing this section, our theory gives a remarkably 
good description of sometimes rather complicated and 
unexpected features of a wide variety of extended ver- 
sions of the Kondo lattice. Moreover, it allows in all cases 
to extract simple physical pictures in order to understand 
the overall trends. The main inaccuracy concerned the 
positions of the /-like Hubbard bands, but we believe 
that this is not a severe deficiency because this is an ex- 
treme high-energy feature. In the cases of extra intra-cell 
repulsions the Lanczos spectra showed some indications 
of a doping induced reconstruction of the electronic struc- 
ture - this could not be reproduced by our theory. Also, 
superconductivity in the doped 'breathing / model' nat- 
urally could not be described. 



V. CONCLUSION 

In summary, we have presented a 'nearly analytical' 
theory for the Kondo lattice. In simplest terms it may 



be viewed as a Fermionic version of linear spin wave the- 
ory, in that we constructed an effective Hamiltonian de- 
scribing the pair creation and propagation of Fermionic 
charge fluctuations on a strong coupling RVB-type vac- 
cuum. This gives us a systematic way to broaden the 
ionization and affinity states of a single cell into bands, 
reminiscent of the cell perturbation theory developed by 
Jefferson and co-workers |Q . For quite a number of dif- 
ferent versions of the Kondo-lattice model, the theoret- 
ical results for the single particle spectral function were 
found to be in overall excellent agreement with exact di- 
agonalization of small clusters, as far as the dispersion 
of energy and spectral weight is concerened. As a par- 
ticularly encouraging fact the rather nontrivial evolution 
of the spectra with electron density in many cases is re- 
produced in detail. This gives us some confidence that 
despite some uncontrolled approximations the theory is 
already very close to the correct picture. Given its ex- 
traordinary simplicity (all results are obtained by an ele- 
mentary Bogoliubov transformation) we believe that the 
theory easily should be amenable to systematic improve- 
ment to describe as yet neglected processes. 
The emerging picture of the Kondo insulator then is 
quite different from that of a conventional band insu- 
lator, and in fact more reminiscent of a superconductor: 
the ground state at half-filling corresponds to an array 
of (overlapping) local singlets, where each Kondo-spin 
forms a bound state with one conduction electron. The 
system has a single particle gap for precisely the same 
reason as in a superconductor: removing a conduction 
electron from some site i and reinserting it at a remote 
site j breaks two f-c pairs, whence there is an increase 
in energy of twice the binding energy of a single pair. 
The basic idea of the present work therefore is not re- 
stricted to the Kondo lattice. Rather, once such an 'RVB- 
vacuum' has been identified, an essentially analogous 
construction can be carried out also for other strongly 
correlated systems. The most obvious example is the 
doped t — J ladder, where the 'rung singlet RVB state' 
may replace the product of single cell states, and for 
which a very similar construction can indeed been carried 
out ||. 

Perhaps the most important approximation we made is 
the neglect of any excited states of two electrons in a 
cell. This means that we have negelected single-cell states 
where one conduction electron and one /-electron couple 
to a triplet. Such states do in fact have a very small 
excitation energy cx V'^/U. As we have already men- 
tioned, this restriction forbids the propagating charge 
fluctuations to 'radiate off' triplet-like spin excitations 
and thus makes their propagation completely coherent. 
The good agreement with the numerics, and also the rel- 
atively 'coherent' nature of the Lanczos spectra them- 
selves (in which the incoherent high-energy continua fa- 
miliar from the t — J oi Hubbard model Bltl are almost 



12 



completely absent), give us some confidence that despite 
their low energy, the neglect of the singlet triplet excita- 
tions is a good approximation as far as the single parti- 
cle properties are concerned. Moreover, it is quite easy 
to incorporate these singlet-triplet excitations into the 
present theory |2^: the three components of a triplet 
state in cell i can be grouped into an SO (3) vector ti. 
Then, we would model the cell i being in the a*^ triplet 
state by the presence of a Bosonic excitation, created 
by In addition to a term describing the on-site en- 
ergies of these triplet we would obtain terms which de- 
scribe the coupling between the Fermionic charge fluc- 
tuations and the Bosonic spin fluctuations. Their form 
can be inferred from rotational invariance: for example 



) (with (T the vector of Pauli matrices) 



describes the hopping of a hole from i — *■ j while creating 
a spin excitation in cell i; the terms (f^ • • o,] r^-j r 
and {al.j.crr,T'aj .,-,) ■ x t^) represent two different ways 
for a spin excitation to exchange its position with a hole. 
The actual prefactors of these terms can be computed in 
an entirely analogous fashion as the hopping integrals for 
the charge fluctuations themselves, and these terms can 
be incorporated into the formalism by using standard 
Green's function techniques. Details will be reported 
elsewhere . Since already the simplest version of the 
theory apparently gives a quite good description of the 
physics, one may expect that such an extension is actu- 
ally 'convergent' and gives a good description of the spin 
dynamics as well. 

As a final remark, we note that the 'vaccum' in the case of 
the simple Kondo lattice was a unique product of singlet 
states, and that the charge fluctuations did have pre- 
cisely the quantum numbers of electrons or holes. For 
systems with a larger unit cell, however, one can envis- 
age situations where this is very different: one example 
would be a 2-fold orbitally degenerate /-level mixing with 
a nondegenerate c-orbital. For one electron/orbital the 
number of electrons /unit cell then would be 3, and the 
(Fermionic) charge fluctuations would correspond to cells 
with an even number of electrons. Thereby Hund's rule 
coupling between the degenerate /-orbitals would favour 
high-spin states, so that one might actually obtain 'high- 
spin quasiparticles'. Form the good success for the simple 
Kondo lattice one might expect that the present formal- 
ism would continue to give a good description even in 
this more complicated situation. 
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